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FOREWORD 


This  report  has  been  prepared  by  Bolt  Beranek  and  Nevnnan  Inc., 
and  describes  the  work  completed  during  the  period  13  May  1962 
through  15  November  1962  under  Contract  No.  AP  33(657) -9118, 
entitled,  "Dynamic  Response  and  Test  Correlation  of  Electronic 
Equipment."  This  work  Is  covered  under  Project  1309»  "Environ¬ 
mental  Simulation  Investigations  for  Flight  Vehicle  Subsystems 
and  Support  Equipment,"  Task  130904,  "Development  of  Simulation 
Facilities  and  Test  Methods  for  Induced  and  Associated  Environ¬ 
ments." 

This  program  was  monitored  by  the  Environmental  Division, 
Directorate  of  Engineering  Test,  Aeronautical- Systems  Division, 
under  the  direction  of  Mr.  Robert  W.  Sevy,  Project  Engineer. 

An  earlier  report  under  Contract  No.  AP  33(6l6)-7743  cover¬ 
ing  the  pei'lod  of  July  1961  to  April  1962  appeared  as  ASD  Technical 
Documentary  Report  ASD-TDR-62-614  entitled,  ''Dynamic  Response  and 
Test  Correlation  of  Electronic  Equipment." 


ABSTRACT 


The  research  has  centered  on  two  main  categories;  the 
replacement  of  acoustic  excitation  by  direct  mechanical  excita¬ 
tion  of  structures  and  the  response  of  substructures  which  are 
tied  to  a  randomly  vibrating  primary  structure.  Studies  are 
Included  on  the  power  transferred  to  a  plate  by  a  shaker,  the 
variation  In  vibratory  response  at  various  positions  and  fre¬ 
quency  bands,  and  the  energy  sharing  between  randomly  excited 
structures  which  have  mechanical  connections.  Flans  for  further 
research  are  also  Included. 


The  publication  of  this  report  does  not  constitute  approval 
by  the  Air  Force  of  the  findings  or  conclusions  contained  herein. 
It  Is  published  only  for  the  exchange  and  stimulation  of  Ideas. 
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iUNDGM  VIBRATION  STUDIES  OF  COUPLED  STRUCTURES 
IN  ELECTRONIC  EQUIPMENTS 


A.  INTRODUCTION 

The  advent  of  rocketiy  telemetry  and  control  requires  elec¬ 
tronic  packages  and  components  which  must  withstand  severe  acoustic 
and  vibrational  environments.  One  step  in  the  design  of  these 
packages  is  the  estimation  of  acoustic  and  vibrational  environments. 
Another  is  concerned  v;ith  vibration  and  sound  reduction  through  the 
design  of  vibration  and  noise  control  treatments.  Finally,  there 
is  a  proof  testing  requirement  which  leads  one  to  seek  methods  of 
establishing  the  expected  package  or  component  environment  in  a 
laboratory  situation. 

These  requirements  all  dictate  a  need  for  better  understanding 
of  the  laws  governing  energy  transfer  between  randomly  vibrating 
systems.  In  the  following  sections  of  this  report  we  shall  describe 
our  efforts  to  improve  this  understanding.  In  Section  B,  we  out¬ 
line  the  general  problem  of  predicting  and  simulating  response  of 
complicated  structures,  ^ie  then  define  certain  subareas  and  tasks 
which  one  must  study  as  part  of  the  more  general  problem.  In 
Section  C  we  describe  our  studies  and  results  thus  far  on  some  of 
these  tasks.  Finally,  in  Section  D,  we  indicate  some  of  the 
additional  tasks  to  be  performed  during  the  remaining  parts  of 
this  contract. 

B.  DESCRIPTION  OF  OVERALL  DIRECTIONS  AND  PROCEDURES 

An  electronic  package  in  a  random  noise  environment  receives 
its  energy  by  two  primary  paths,  direct  acoustic  excitation  and 
mechanical  vibration  through  its  attacliment  points.  We  will 
normally  assume  that  the  power  from  these  two  sources  is  additive. 
Correlation  effects  which  would  cause  departures  from  additivity, 
are  possible  since  the  surrounding  structure  generates  sound 
about  the  package.  The  vibrations  of  the  base  structure  may 
therefore  be  correlated  with  the  sound  field.  The  difference  in 
propagation  times  for  sound  and  mechanical  energy  and  the  differ¬ 
ences  in  excitation  efficiencies  may  be  expected,  however,  to 
minimize  the  effect  of  any  correlated  behavior. 
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If  one  ccxnputed  the  power  absorbed  by  the  package  from  the 
sound  field  and  ^•hrouj^  its  mounting  points,  then  a  proof  test 
of  the  design  can  be  made  by  the  direct  application  of  this  power. 

To  apply  this  power  with  a  sound  field  however  requires  very  expen¬ 
sive  sound  sources.  It  has  been  suggested  that  one  might  be  able 
to  substitute  direct  mechanical  excitation  instead. 

The  power  fed  into  a  panel -frame  structure  by  a  random  noise 
sound  field  has  been  studied  at  BBN  both  theoretically  and  experi¬ 
mentally  (Refs.  1-3).  The  approach  is  statistical  in  its  description 
of  the  random  noise  process,  and  it  is  also  statistical  in  its 
description  of  the  structure.  Rather  than  analyze  a  structure  with 
precise  dimensions,  mounting  conditions,  damping  treatments,  etc.., 
one  describes  the  structure  by  an  average  modal  density,  total 
area  and  edge  length  of  panels,  loss  factor,  etc.  This  allows  one 
to  apply  the  results  to  a  wider  class  of  structures,  but  one  must 
also  be  concerned  with  how  variations  in  response  will  arise  from 
one  structure  to  another  due  to  the  variations  in  these  statistical 
measures.  Ihese  same  considerations  will  arise  when  we  discuss 
how  energy  is  transferred  between  structural  attachments. 

In  ttils  present  program,  we  have  asked  the  following  specific 
questions: 

a)  What  are  the  laws  of  energy  transfer  between  connected 
structures?  The  answer  to  this  tells  us  how  much 
energy  will  get  to  a  package  through  its  mountings  and 
further,  the  amount  of  energy  which  travels  to  internal 
panels  of  the  package. 

b)  Since  our  answers  to  (a)  are  essentially  average  values 
of  statistical  results,  what  variations  may  we  expect 
as  we  go  from  one  structure  to  another,  or  from  one 
part  of  the  frequency  spectrum  to  another,  or  from  one 
position  on  the  structure  to  euiother?  A  knowledge  of  the 
variations  will  allow  us  to  set  confidence  levels  on  any 
range  of  estimates  we  may  wish  to  consider. 

c)  If  we  wish  to  supply  mechanical  power  to  the  package, 
not  through  its  normal  environmental  condition,  but  by 
diz*eot  application  of  high  frequency  shakers,  can  we 
develop  reliable  estimates  of  input  power?  Also,  can 

we  develop  techniques  for  measuring  the  input  mechanical 
power  directly? 

d)  Finally,  can  we  use  these  theoretical  results  to  actually 
predict  response  of  a  simple  package -sub -panel  system  in 
a  noise  and  vibration  environment?  If  the  answer  to  this 
is  "ycBf"  then  our  confidence  in  the  theoretical -experi¬ 
mental  procedures  referred  to  above  is  augmented. 
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C.  DESCRIPTION  OP  SPECIFIC  STUDIES  AND  RESULTS 

We  will  now  describe  specific  studies  which  we  have  launched 
to  answer  the  questions  posed  in  the  preceding  section.  These 
studies  are  not  complete^  but  we  feel  that  a  good  beginning  has 
been  made.  The  studies  themselves  appear  as  a  series  of  Appendices 
to  the  report.  Here  we  shall  merely  review  their  motivation  and 
discuss  the  main  results  very  briefly. 

I.  Power  Flow  Into  Structures  Excited  Mechanically 

A  shaker  may  be  mounted  to  a  structure  by  a  stud  attachment. 
Normally  this  stud  is  sufficiently  small  compared  to  a  wavelength 
of  bending  waves  and  the  structural  dimensions  that  it  may  be 
considered  a  point  force.  The  power  this  source  feeds  into  the 
structure,  therefore,  depends  on  the  point  impedance  at  the  posi¬ 
tion  of  application.  This  impedance  is  usually  a  highly  fluctuating 
function  of  frequency  for  finite  stz*uctures,  but  if  one  is  con¬ 
cerned  with  average  response  over  a  band  of  frequencies,  then  the 
infinite  system  input  impedance  may  theoretically  be  used. 

The  study  of  a  shaker  driven  plate  in  Appendix  1  is  an  attempt 
to  test  the  assumption  that  the  average  power  fed  into  the  plate 
in  third  octave  bands  can  be  calculated  on  the  basis  of  the  knov/n 
force  from  the  shaker  and  the  input  resistance  of  the  infinite 
plate.  The  early  conclusion  seems  to  be  that  this  is  possible  if 
one  will  accept  an  uncertainty  of  the  order  of  3  db. 

II.  Response  Variations  in  Multimodal  Structures 

Most  of  the  studies  of  sound  to  vibration  and  vibration  to 
vibration  Interaction  generate  estimates  of  average  energy  levels 
or  power  flow  in  frequency  bands.  These  estimates  are  more  nearly 
exact  as  one  averages  over  more  modes  of  response,  positions  of 
observation  and  types  of  sti*uctures.  Some  concern  has  arisen 
about  the  confidence  one  can  have  in  these  average  values  if  one 
is  dealing  with  only  a  few  modes,  or  observation  points,  or  struc¬ 
tures.  To  generate  directly  procedures  for  evaluating  confidence 
levels  for  the  estimates,  we  have  initiated  a  set  of  studies  on  the 
variations  in  response  for  various  cases. 

In  Appendix  II  these  variation  studies  are  reported.  We  have 
studied  the  variations  in  response  at  one  position  when  the  struc¬ 
ture  is  point  excited  at  a  second  position  (Sections  B  and  C)  and 
the  response  of  a  simple  oscillator  when  it  is  attached  to  the 
plate  (Section  D) .  The  first  study  is  appropriate  to  expected 
variations  in  estimating  the  vibration  of  the  plate  itself  while 
the  second  is  concerned  with  variations  in  substructure  response  — 
in  particular  in  regions  of  low  modal  density. 
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Ill*  Energy  Sharing  Between  Connected  Structures 

The  power  flow  between  structures  and  their  vibration 
levels  are  studied  by  considering  each  structure  as  a  collection 
of  modes.  These  modes  are  coupled  by  the  connection  of  the 
structures.  The  evaluation  of  this  coupling  requires  detailed 
calculation  of  Input  Impedances  of  various  types  for  Infinite 
systems  Just  as  the  sound-structui*e  problem  requires  computation 
of  radiation  resistance  for  various  configurations.  Our  work  on 
connected  structuz*es,  therefore,  has  Involved  the  computation  of 
these  Impedances  for  cantilevered  beams  attached  to  plates  and 
for  plates  Joined  together  either  by  continuous  welds  or  studs. 

It  has  also  involved  experiments  on  the  svstems  analyzed  theoreti¬ 
cally  to  establish  whether  In  fact  one  can  get  agreement  with 
theoretical  estimates. 

The  details  of  the  work  In  this  area  are  reported  In 
Appendix  III.  In  general,  our  results  for  vibration  levels  In 
the  cantllevez%d  beam  and  the  welded  plate  are  quite  encouraging. 

D.  CONCLUSIONS 

Power  flow  Into  stroctures  excited  mechanically.  In  this 
area  we  are  developing  techniques  for  measuring  the  force  velocity 
product  so  that  Input  power  may  be  measured  directly  for  pure 
tones  or  bands  of  noise.  This  means  that  experimentally  we  do 
not  have  to  rely  on  computation  of  Input  Impedances  or  structures 
with  the  attendant  uncertainty  In  parameter  values. 

Variance  analysis  of  multimodal  structures.  We  now  feel 
that  variance  and  confidence  analysis  has  been  sufficiently  well 
established  as  a  theoretical  tool.  We  do  not,  therefore,  plan 
further  theoretical  studies  In  developing  the  method.  It  will 
Instead  become  a  part  of  the  power  flow  and  energy  sharing  studies 
and  be  used  and  tested  In  the  experimental  aspects  of  that  woz4c. 

Energy  sharing  between  connected  structures.  We  should  like 
to  consider  at  least  two  systems  which  have  not  been  studied 
experimentally  at  this  time.  One  Is  a  study  of  plates  connected 
by  stud  mountings,  a  type  of  connection  which  Is  very  common  In 
real  packages.  The  second  Is  the  response  of  a  single  oscillator 
attached  to  a  plate.  This  is  essentially  a  test  of  some  of  the 
variance  studies  imported  In  Appendix  II.  In  addition  to  these, 
some  additional  theoretical  studies  of  plate -beam  decoupling  by 
notching  of  the  beam  would  be  In  order. 
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Environmental  estimates  and  Bimulatlon  for  package-panel 
structure.  We  plan  to  take  a  box  with  an  internal  panel  and  place 
It  in  a  sound  field,  measuring  the  resulting  vibrational  levels. 
With  a  knowledge  of  the  damping,  we  expect  to  correlate  these 
levels  with  the  radiation  resistance  of  the  structure.  This  part 
of  the  experiment  Is  essentially  a  repeat  of  am  earlier  study 
(Ref.  3).  Then  we  will  apply  some  Input  power  with  a  shaker  and 
again  measure  vibration  levels,  with  special  attention  to  seeing 
whether  the  acoustic  response  can  be  simulated  In  this  way. 
Comparisons  of  internal  panel  response  with  package  wall  response 
should  afford  a  reasonable  test  of  our  energy  sharing  calculations. 


5 


NOTES  AND  REFERENCES 


1.  R.  H.  loTon  and  G.  Maldanlk,  "Power  flow  between  linearly 
coupled  oscillators,"  J.  Acoust.  Soc.  Am.  3^»  5j  623, 

May,  1962. 

2.  G.  Maldanlk,  "Response  of  ribbed  panels  to  reverberant 
acoustic  fields,"  j.  Acoust.  Soc.  Am.  6,  809,  June,  1962 

3.  D.  U.  Nolseux,  J.  J.  Coles,  N.  Doelllng,  R.  H.  Lyon, 
"Response  of  electronics  to  Intense  sound  fields,"  Report  No 
871 J  Bolt  Beranek  and  Newman  Inc.,  Cambridge,  Massachusetts, 
July,  1961,  Section  III. 


6 


APPENDIX  I 


POWER  FLOW  INTO  STRUCTURES  EXCITED  MECHANICALLY 


A.  INTRODUCTION 

In  this  appendix  we  consider  the  injection  of  mechanical 
power  into  a  structure  by  an  electromechanical  shaker.  The  first 
step  is  the  determination  of  the  coupling  between  the  mechanical 
shaker  (or  a  system  of  shakers)  and  the  structure  we  wish  to  test. 
We  propose  to  substitute  this  mechanical  excitation  for  the 
acoustic  power  supplied  to  the  structure  by  a  sound  field. 

The  power  absorbed  from  a  reverberant  sound  field  of  mean 
square  pressure  ^  is  given  by  (Ref.  1) 


where 

Mg  is  the  mass  of  the  structure, 

a)  is  the  center  (radian)  frequency  of  excitation, 

c  is  the  speed  of  sound  in  air, 

Ap  is  the  irradiated  area  of  the  structure, 

ng(u>)  is  the  average  modal  density  of  the  structure, 

a  d  radiation  efficiency  of  the  structure 

averaged  over  Ap  and  the  frequency  band  of 
the  excitation. 

It  is  this  power  that  v/e  want  to  supply  by  the  shaker  system. 
The  power  which  a  shaker  can  feed  to  the  structure  will  depend  on 
the  real  part  of  the  mechanical  impedance  presented  to  it.  We  may 
define  the  mechanical  impedance 


^in  =  «in  -  ^  *ln 


(1.2) 


The  power  Input  In  terms  of  the  mean  square  acceleration  at  . 
the  shaker  ag  (in  the  absence  of  reverberation)  is 
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(1.3) 


where  we  have  assumed  that  Is  a  fairly  smooth  function  of 
frequency.  In  the  initial  stages  of  the  computations  we  shall 
use  an  Impedance  which  is  averaged  over  several  modes  of  vibra¬ 
tion.  If  we  equate  Equation  (l.l)  with  Equation  (1.3)  we  obtain 
for  the  meaui  square  pressure 


a2  ^inj^s 

^  2^^"s(“)c%Orad 


(1.4) 


We  shall  be  studying  the  mechanical -acoustical  equivalence  in  the 
case  of  panels.  For  a  panel  of  thickness^  h,  material  density, 
Pp,  and  longitudinal  velocity, c^,  the  modal  density  is  given  by 

(Rer.  1) 


ng(ai) 


\/3  As 

Sirhc 

Ai 


(1.5) 


where  Ag  is  the  panel  surface  area  (one  side). 

The  expression  for  the  impedance  for  an  infinite  panel  is 


For  a  finite  panel  the  same  equation  is  obtained  vAien  the  input 
impedance  is  averaged  over  several  modes.  Our  assumption  is  l^t 
in  a  narrow  frequency  band  (e.g.,  third  octave)  several  modes 
of  the  pwel  will  be  excited. 
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The  radiation  efficiency  for  a  baffled*  panel  is  (Ref.  2) 


°ra.d 


IT  h  c  ■ 
r  * 


(1.7) 


where  Pj.  is  the  perimeter  of  the  panel  and  fp  is  the  acoustic 
coincidence  frequency.  The  expression  for  is  valid  if  the 

sound  field  and  the  vibrational  field  of  the  structure  are 
reverberant  and  the  dimensions  of  the  panel  are  greater  than  one- 
third  of  an  acoustic  viavelength.  For  other  frequency  conditions, 
the  exprecslon  is  somewhat  more  complicated,  (Ref.  2).  Substitu¬ 
ting  Equations  (1.5)j  (1.6),  and  (1.7)  in  Equation  (1.4),  we 
obtain 


PpJ^ 

VT 


(1.8) 


For  a  Goodmans  shaker  (Goodmans  Industries,  Ltd.,  Model  390A) 
the  limiting  rms  acceleration  v/hen  unloaded  is  about  100  g*s.  If 
the  shaker  is  attached  to  a  point  of  a  panel  vjhose  point  impedance 
is  the  mean  square  acceleration  at  the  point  is 


(1.9) 


v/here  can^  is  the  mass  reactance  of  the  moving  coil  assembly.  This 
is  the  mean  square  acceleration  which  we  use  in  defining  Pin 
Equation  (l^S).  For  a  panel  with  h  =  3.16  x  lO'-l  cm,  pp  =  2.7 
Cj  s  5  lo5  civ^sec,  auid  Pr  =  6l0  cm,  the  maximum  available  equiva¬ 
lent  sound  pressure  level  in  db  is  about 


In  this  section  we  are  interested  in  determining  the  power 
fed  into  a  stinicture  by  mechanical  excitation.  For  the 
specific  case  of  a  panel  excited  by  mechanical  means,  whether 
the  panel  is  baffled  or  not  is  not  important,  as  long  as  the 
radiation  damping  is  smaller  than  the  mechanical  damping. 
However  when  one  wants  to  study  experimentally  the  equivalence 
betv/een  the  acoustical  and  mechanical  power  supplied  to  the 
structure,  whether  the  panel  is  baffled  or  unbaffled  becomes 
essential.  For  an  unbaffled  panel  the  radiation  efficiency 
^rad  lower  than  the  one  given  in  Eq.  (1*7). 
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142  +  5  logj^o  ^  ^°®10 


(1.10) 


(«c)^  ■.  »L 

For  the  panel  under  consideration  f p  <2:  ^  (aluminum  panel) 
and  -  3.12  x  10^  gny^sec.  The  mass  of  the  moving  coll  assembly 

of  the  Goodmans  shaker  (Model  390A)  Is  about  72  gm.  Thus,  the  maxi¬ 
mum  equivalent  sound  pressure  capability  of  this  Goodmans  shaker 
with  respect  to  the  given  panel  Is  at  100  cps  about  138  db  and 
at  1  kc  about  l48  db. 

VIhen  an  accelerometer  Is  placed  at  the  point  to  which  the 
shaker  Is  connected  the  acceleration  measured  at  this  point  may 
be  influenced  by  the  reverberant  field  of  the  panel.  The  exact 
expression  Is 


>1 


Q  r 


'in 


'In 


2  2 
o)  ra„ 
c 


(1.11) 


where  a  Is  the  mean  square  acceleration  Indicated  by  the  acceler¬ 
ometer  Snd  a?  Is  the  mean  square  acceleration  of  the  reverberant 
vibrational  Afield. 


The  power  supplied  by  the  shaker  must  be  equal  to  the  power 
dissipated  In  the  system.  This  equality  enables  us  to  test  our 
formalism  experimentally.  The  damping  In  the  panel -shaker  system 
may  be  expressed  In  terms  of  a  loss  factor 


The  dissipated  power  Is 

*’dlss  ^”s  ’V'“>)  (1.13) 


where  Tp  Is  the  i^everberatlon  time  of  the  system.  Equating 
Equation  (1.3)  to  Equation  (1.13)f  and  making  use  cf  Equation 
(1.11),  so  that  the  meeui  square  accelerations  Involved  are 
measurable  quantities,  we  obtain 
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+ 


(1.1^) 


s  r 


Mglp 


'in 


—K  ~K 

Thus,  measurements  of  will  give  us  the  experimental  value 

of  the  bracketed  quantity  In  Equation  (l.l4).  On  the  other  hand, 
measurements  of  q  Tby  reverberation  time  measurements  and  the  use 
of  Equation  (1.12) j,  the  use  of  the  given  value  of  me  >  72  gm  In 
our  case)  and  the  theoretical  value  of  Rin  &&  given  by  Equation 

il.6)  will  enable  us  to  calculate  the  bracketed  factor  In  Equation 
1.14)  by  an  Independent  method.  Comparison  between  these  two 
determinations  should  give  us  a  cross  check  for  our  theory.  It 
may  be  In  order  to  summarize  at  this  point  the  fundamental  assump¬ 
tions  that  we  made  and  that  therefore  must  be  satisfied  for  fair 
comparison  to  be  possible. 


a)  The  vibrational  field  of  the  panel  is  reverberant 
In  the  sense  that  all  or  most  of  the  modes  In  a 
given  frequency  band  (e.g.,  third  octave)  Are  parti¬ 
cipating  In  the  vibration  and  that  the  vibrational 
field  Is  diffuse.  An  experimental  test  for  the 
reverberant  field  Is  that  the  vibrational  level  of 
the  panel  when  excited  by  noise  filtered  thz^iugh  a 
given  frequency  band  will  be  essentially  Independent 
of  position. 

b)  The  point  Impedance  of  the  panel  approximates  that 

of  the  point  Impedance  of  an  infinite  pemel.  Theoreti¬ 
cally,  If  (a)  Is  satisfied,  so  is  (b). 

B.  EXPERIMENTAL  TEST  OP  EQUATION  (l.l4) 

An  experimental  arrangement  has  been  designed  and  constructed 
to  test  Equation  (I.l4)  which  Is  derived  In  Section  A.  The  differ¬ 
ence  In  acceleration  levels  at  the  driving  point  and  a  selected 
number  of  points  on  the  panel  were  measured.  The  reverberation 
time  of  the  panel -shaker  system  was  also  determined  experimentally 
and  the  bracketed  expression  In  Equation  (l.l4)  was  calculated  In 
accordance  with  the  procedure  Indicated  at  the  end  of  Section  A. 

The  latter  calculations  vere  then  compared  with  the  measured 

-vrv 

values  of  the  ratio  results  for  three  different  cases, 

undaunped,  lightly  damped  and  heavily  damped  panels,  are  given  In 
graphical  form  In  Figures  1.1,  1.2,  and  1.3> 
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1.  Experimental  Apparatus 


An  aliimlnum  panel  3'  x  4*  x  0.129"  was  suspended  from  the 
celling  on  nylon  cords  from  two  comers.  Otherwise  It  was 
unrestrained.  A  vibration  shaker  was  arranged  to  excite  the 
panel  at  a  point.  The  location  of  the  driving  point  was  chosen 
away  from  the  edge  of  the  panel.  The  attaching  hazKlware  was 
devised  to  provide  only  treunsverse  motion  and  minimize  applied 
torque  and  damping  at  the  driving  point. 

'Oie  excitation  was  in  third  octave  bands  of  random  noise. 

The  output  of  tftie  accelerometer  attached  at  the  driving  point 
was  compared  successively  with  the  output  of  the  accelerometer 
attached  at  several  points  on  the  paunel  surface  (see  Figures  1.1, 
1.2,  and  1.3  for  the  exact  locations  which  are  Indicated  by 
circled  numbers).  The  output  of  the  accelerometers  was  recorded 
on  a  graphic  level  recorder.  The  flow  diagram  of  the  experimental 
arrangement  Is  shown  In  Figure  1.4. 

The  reverberation  time  was  determined  using  a  graphic  level 
mcorder  and  the  reverberation -time  apparatus  described  In 
reference  4  of  Appendix  III.  A  pulse  technique  was  used  for  the 
reverberation  time  measurements;  the  panel  was  tapped  with  a 
metallic  hammer  at  various  positions  on  the  panel  and  Its  decay 
was  measured  by  measuring  the  filtered  (third  octave  band)  output 
of  an  accelerometer.  The  flow  diagram  for  the  experimental 
arrangensot  Is  shown  In  Figure  1.5. 

2.  Discussion 

Althoutfi  there  Is  a  general  correlation  between  the  results 
of  ^e  two  methods  of  determining  the  factor  of  proportionality 

relating  a^  to  a^.  In  detail  the  agreement  Is  not  completely 

satisfactory.  Further  Improvement  In  the  theoretical  approach 
may  therefore  be  necessary.  It  should  be  pointed  out  that  for 
the  undamped  and  lightly  damped  case  the  vibrational  field  could 
be  assumed  reverberant.  The  measurements  taken  at  the  ttu?ee 
positions  on  the  pauiel  did  not  vary  by  more  than  2  db.  In  the 
case  of  ttw  heavily  damped  peuiel  the  variation  In  level  between 
the  four  posltlcwis,  where  measurements  were  taken.  In  several 
Instances  exceeded  6  db.  These  large  differences  In  level  were 
not  confined  to  any  particular  frequency  range.  Moreover,  the 
measurements  of  reverberation  time  for  the  heavily  damped  panel 
at  frequencies  above  3  kc  are  not  very  reliable;  they  were  taken 
at  the  limit  of  the  range  of  the  Instrument  used.  The  reverbera¬ 
tion  times  above  3  ke  were  of  the  order  of  10-2  sec.  In  Mils 
ease,  rfe  note  here  that  some  of  the  theoretical  results  from 
Section  B  In  Appendix  11  anticipate  the  hlghoz*  fluctuations  In 
response  for  ttie  very  heavily  damped  panel. 
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FI6.I.I  ACCELERATION  LEVEL  DIFFERENCE  BETWEEN  DRIVING  POINT 
AND  REVERBERANT  FIELD  FOR  3'x  4'  PANEL 
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FIG.  1.2  ACCELERATION  LEVEL  DIFFERENCE  BETWEEN  DRIVING  POINT 
AND  REVERBERANT  FIELD  FOR  3'x4'  PANEL 
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FI6.L3  ACCELERATION  LEVEL  DIFFERENCE  BETWEEN  DRIVING  POINT 
AND  REVERBERANT  FIELD  FOR  3'x4'  PANEL 


FIG.  1.4  PANEL  EXCITATION  SYSTEM 
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APPENDIX  II 


RESPONSE  VARIATIONS  IN  MULTIMODAL  STRUCTURES 


A.  INTRODUCTION 

In  the  past  few  years,  the  response  of  complex  structures 
made  up  of  panels,  supporting  members,  struts,  etc.,  to  random 
noise  fields  has  received  a  good  deal  of  attention  at  Bolt 
Beranek  and  Newman  Inc.,  and  other  research  centers.  At  BBN  we 
have  described  not  only  the  temporal  response  In  a  statistical 
way,  but  the  structure  itself  as  well  (Ref.  1).  Thus,  a  struc- 
ttire  is  defined  by  its  distribution  of  modal  resonaunces  in  frequency, 
average  lengths  of  beams,  areas  of  panels,  etq.  The  results  of 
this  procedure  have  been  encouraging  in  their  ability  to  predict 
certain  average  characteristics  of  response  lAien  many  modes  of 
vibration  are  involved. 

For  small  structures  in  the  lower  frequency  ranges,  one  may 
not  always  have  a  ve3?y  large  number  of  modes  contributing  to  ttxe 
response.  One  effect  of  this  will  be  to  cause  fluctuations  about 
ttie  average  value  in  the  observed  response  from  one  frequency 
band  to  the  next  (Ref.  2),  which  we  shall  refer  to  as  frequency 
variation.  Also,  there  will  be  variation  in  response  as  one 
moves  from  one  measuring  point  to  another  on  the  structure.  This 
spatial  variation  will  be  Important  as  an  additional  source  of 
variation  in  the  comparison  of  theoretically  computed  average 
values  and  experimentally  measurad  values.  Also,  one  should  note 
that  we  are  not  concerned  here  with  the  temporal  fluctuations  of 
the  signal  as  a  random  process  in  time.  We  assume  that  sufficient 
averaging  time  is  used  for  the  measurement  so  that  variations  due 
to  UU.S  effect  are  negligible. 

It  Is  ambiguous  to  speak  merely  of  "the  variation"  since  Uie 
response  variable  we  have  In  mind  has  not  been  specified.  Analyt¬ 
ically,  ^e  variation  in  (time)  average  response  (velocity, 
acceleration,  strain,  etc.)  is  convenient  to  deal  with.  Experi¬ 
mentally,  response  is  read  out  and  plotted  In  logarltlxnic  units, 
and  the  variation  visible  on  the  experimental  plot  is  in  decibels, 
not  so  many  g's  (for  example).  Unfortunately,  a  great  deal  must 
be  known  about  the  response  in  oz^ler  to  get  its  logarithmic  variation. 
If,  for  example,  we  know  the  standard  deviation  of  the  acceleration, 
we  cannot  go  from  this  value  to  a  standard  deviation  for  the 
response  in  decibels.  Indeed,  the  entire  probability  density  for 
the  response  would  be  required  to  evaluate  ttie  logarithmic  variation. 
In  some  cases  we  shall  find  it  possible  to  make  this  computation, 
in  others  we  will  have  to  be  satisfied  with  the  variation  in  the 
dynamical  variable  itself. 
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B.  1»E  VARIATION  OF  RMS  RESPONSE  IN  BANDS  OF  NOISE 
FOR  MULTIMODAL  SYSTEMS 

We  wish  to  find  the  variation  In  mean  square  response  (veloo- 
ItYi  say)  for  a  stinicture,  which  has  a  pure  tone  mean  square 
response*  given  by  W((o),  when  that  structure  Is  excited  by  a  band 
of  noise  of  width  A  centered  about  the  frequency  o.  The  response 
for  band  excitation  of  bandwidth  A  and  spectral  density  ^  Is 


V(“o) 


1 


W((u)da). 


(2.1) 


Let  the  average  response  (TCo)),  be  the  response  W((u)  averaged 
over  similar  structures  or  over  different  positions  on  the  same 
structure.  The  b^d  response,  to  excitation  of  bandwidth  A  and 
spectral  density  ^  ,  of  the  average  response  ^a>)  Is 


^^+V2 

J%-V2 


(2.2) 


provided  7((o)  varies  only  slightly  over  the  bandwidth  of  excitation 
due  to  changes  In  modal  density  or  damping. 

One  measure  of  the  variation  In  response  for  band  excitation 
Is  the  standard  deviation  of  v  which  Is  found  from  the  varlamce 


’2  jJJn 


(2.3) 


dcUn  P 
2  w 


(«>l-<»2) 


"z^sponse"  In  this  section  means  "mean  square  response 
In  time"  in  all  cases. 
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where 


o 


2 

w 


Is  the  variance  of  the  sinusoidal  3?esponse, 
and 


o 


w 


(2.4) 


In  the  case  of  multimodal  systems  with  a  hi^  degree  of  modal 
overlap  (a  situation  appropriate  to  rooms  but  not  as  common  in 
complex  structures)  Schroeder  has  shown  that  the  correlation 
function  for  pure  tone  response  (Ref.  3)  is  of  the  form 


«w(0 


4a^ 

^  +  4a^ 


(2.5) 


where  a  is  the  modal  decay  constant,  related  to  the  structural 
loss  factor  R  by 


ij  =  20/(0  .  (2.6) 

Also,  Schroeder  has  shown  (Ref.  4)  that  the  distribution  of 
W  is  given  by  the  exponential  density 


,  -W/nV, 

M  e  Mw  >  o) 


(2.7) 


which  has  a  variance 


2  2 
®w  =  ^w 


(2.8) 


If  we  place  Equation  (2.5)  in  Equation  (2.3)  and  transform 
to  variables 
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s  -  0)^  -  (ttg 

TJ  -  0)^  +  (Dg  \ 

d4  dT)  -  2  dto^  diOg 

/ 

then  the  Integration  becomes 


(2.9) 


a 


2 

V 


N  ®w 


- ^  in  (1  +  ir^^) 

ffTJ 


N  ®w 


f(N) 


(2.10) 


where  N  •>  n/^aat  Is  the  ratio  of  excitation  bandwidth  to  the  effec¬ 
tive  bandwidth  of  the  resonant  modes.  As  N— — )<»,  the  tem  In  the 
brackets  which  Is  plotted  In  Figure  2.1  goes  to  unity  and  one  has 
a  result  reminiscent  of  a  central  limit  type  of  convergence  of  V 
on  Its  mean  value. 


We  note  that  the  ratio  between  band  variance  Oy  and  pure  tone 
variance  ow  depends  only  on  the  ratio  of  exciting  bandwidth  to 
modal  bandwidth  only;  there  Is  no  dependence  on  modal  spacing 
because  of  the  model  used.  This  may  be  because  we  have  assumed  a 
high  degree  of  modal  overlap  In  using  Schroeder's  formulas,  which 
means  that  the  modal  bandwidth  2ra  is  large  compared  to  the  average 
spacing  between  modes.  In  Section  C  of  this  Appendix,  we  extend 
these  results  to  the  region  of  low  modal  density  which  Is  a  more 
appropriate  z^glme  for  structures. 
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We  note  from  the  result  Equation  (2.10)  that  the  fluctuation 
Increases  as  the  damping  a  increases.  Our  Intuitive  feeling,  how¬ 
ever,  Is  that  the  response  should  "smooth  out"  as  damping  is  added, 
nils  Is  true  In  the  low  mode  overlap  region  to  be  studied.  For 
high  modal  overlap  however,  the  rapidity  of  fluotuatlcm  In  tiie 
pure  tone  response  decreases  with  Increased  damping,  while  the 
variation  In  response  remains  more  or  less  constant.  Thus  If  we 
average  this  response  over  a  fixed  band  of  frequencies,  we  average 
over  fewer  variations  when  the  damping  Is  higher,  and  we  may  expect 
to  see  greater  variation  In  the  average. 

The  variance  of  band  response  can  be  expressed  in  terms  of  Its 

mean 

f(N)  =Hr(N) 


-  I  f(N)  mj  , 

or. 


where  we  have  used  Equation  (2.2) • 

Thus,  for  the  average  Vp  of  a  set  of  measured  responses. 


1.  Example 

_2 

A  large  panel  has  a  loss  factor  q  =  10  and  an  average  modal 
separation  of  4  cps.  A  resonant  filter  Is  equivalent  to  a 
rectangular  filter  of  bandwidth  vqf.  At  1  kc,  this  effective 
modal  bandwidth  Is  31.4  cps.  Thus,  on  the  average  there  Is  an 
overlap  of  7  or  8  modes  which  Is  probably  sufficient  to  meet 
Schroeder's  criterion.  We  can  ask  what  the  standard  deviation 
would  be  In  the  third  octave  band  centered  on  1  kc.  In  this 
case,  the  bandwidth  of  excitation  Is  260  cps  amd 

N  =  “  8-28  • 
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Prom  Figure  2.1,  f(8.28^  ■  0.9  end 


a 


Thus,  ttie  etandaz^  deviation  In  reeponee  will  be  about  3  db 
from  the  average.  ^  repeating  the  measurement  at  several  points 
on  ttie  structure,  tMs  uncertainty  may  of  course  be  reduced.  If 
we  take  t)  -  2  x  10~3,  then  vT)f  -  6.3  ops,  and  N  -  4l.4.  Thus 


which  Is  only  about  1  db  variation  about  the  mean.  Of  course,  we 
are  now  probably  violating  the  condition  of  strong  modal  overlap, 
but  the  point  can  be  made  that  damping  Increases  the  fluctuation 
In  this  region.  This  observation  appears  to  be  supported  by  some 
of  our  preliminary  experiments  (see  Appendix  1). 

C.  THE  VARlATlCm  OF  RMS  RESPONSE  IN  MULTIMODAL  STRUCTURES 

VIHEN  MODAL  OVERLAP  IS  ^LL 

In  Section  B  we  showed  that  the  band  response  of  the  structux>e 
was  related  to  the  pure  tone  response  In  a  very  simple  mauiner.  In 
particular,  the  variance  of  band  and  pure  tone  response  were 
related  as  Indicated  In  Eqiaatlon  (2.3).  This  relation  Is  unaffected 
by  any  assumptions  of  modal  density,  deunplng,  etc.  We  may  therefore 
use  It  as  a  starting  point  In  our  present  analysis. 

Again  calling  W((o)  the  mean  square  pure  tone  i^sponse  we 
require  Its  autocorrelation  function  0w(s).  This  Is  found  by 
first  considering  the  response  aunplltude 


where  bm  Is  an  excitation  coefficient  and  gm(o>)  Is  a  harmonic 
oscillator  (modal)  response  function. 


24 


(3.2) 


Here,  (u^  Is  the  modal  resonance  frequency  and  Om  Is  related 
to  the  modal  loss  factor  by  ■  \* 

If  W  Is  the  metm  square  pure  tone  response,  then 

W(<o)  -  g-  ly(aj)  1  .  (3.3) 


Ae  shall  assume  that  all  modes  are  fairly  lightly  damped 
so  that  strong  resonances  occur  at  a>  »  oc^.  The  expression  for 
W  may  be  written 


m^n 

(3.4) 


m 


where  we  have  assumed  that  ginSn**”^,  when  cum  /  (cu  for  any  cu. 

This  Is  essentially  the  assumption  of  very  small  mo^l  overlap. 

The  expression  for  W  Is  a  superposition  of  "events" lg„,^  of 
12 

magnitude  g-  b^  spaced  along  the  frequency  axis.  For  a  fiat  panel, 

the  average  modal  density  Is  constant  (Ref.  3),  &nd  If  we  assume 
a  Poisson  placement  of  modes  along  the  axis,  we  can  directly 
compute  some  needed  statistical  properties  of  W.  We  shall  also 
assume  that  the  damping  Is  the  same  for  all  modes,  ■  a  >  const. 

Rice  has  shown  that  the  correlation  function  of  a  process 
like  W  has  a  correlation  function  given  by  (Ref.  6) 


'^„(4)  ng(a3)  g^-#g^  ,  (3.5) 
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where  Is  the  symbol  for  the  convolution  operation  defined 

by 


g 


2  2 

l8(“)  I  ls(uH-0  1 


(3.6) 


Since  |g|^  has  a  spectrum  given  by^ 

Sg(T)  -  e  , 

the  spectrum  of  the  convolution  of  |g|^  Is  (Ref.  7) 

0(T)  - 


nie  Inverse  of  this  spectrum  Is  the  convolution. 


Isl^*  Ul=  - 


Placing  this  In  Equation  (3.5)  and  requiring  that  - >1 

as  i - >0  gives 


♦w(0 


(3.7) 


The  result  Equation  (3.7)  la  the  same  as  that  which  Schroeder 
obtained  for  the  case  of  high  mode  density  given  In  Equation  (2.5). 
This  strongly  suggests  that  the  correlation  function  for  4  Is  not 
dependent  on  modal  spacing.  Since  the  correlation  function  Is 
given  by  Equation  (3.7)«  z*esult  of  Equation  (2.10)  for  the 
varlaunce  of  band  response  Is  unchanged.  The  expression  for  m„  and 
0  will  be  altered  however. 
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The  mean  value  of  w  Is  (Ref.  8) 


ST 


The  variance  of  W  Is  given  by  (Ref.  9) 
®w  “  ?  "s^"^  ^  f  dm  |g(tB)  1^ 


(3.8) 


(3.9) 


The  excitation  coefficients  b  depend  on  the  type  of  source 
used  to  drive  the  structure.  If  we  use  a  shaker  attached  at  a 
point  x®  they  are  of  the  form 

b  ~  B  sin  mj^vXj^/ij^  sin  mg^Xg/^^g  sin  sin  mgVX^ig  (3.10) 


We  may  form  the  moments  of  b  by  averaging  over  x,  x  : 

laT 


Thus, 


~T  4 
b  - 


i]  - 


™  m  r, 

™w  “  "V  “  sr  "S  ^ 


"S  -  (l| 


■I  "i? 

5iran, 


(3.11) 


Using  Equation  (3.11)  In  Equation  (2.10),  we  have  finally. 
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(3.12) 


We  note  from  Equation  (3.12)  that  as  the  ntamber  of  modes  naA 
embraced  by  the  exciting  bandwidth  Is  Incz^ased,  the  variation 
In  response  Is  diminished.  This  agrees  with  our  Intuition.  Also, 
as  the  damping  Is  diminished,  N  Increases  and  f(N)  Increases, 
rapidly  for  small  N  and  more  slowly  as  N  rises  above  unity  ( see 
Figure  2.1).  Thus,  the  variation  In  response  Increases  as  damping 
Is  reduced,  ac^ln  In  agreement  with  our  Intuitive  notions.  The 
behavior  In  the  low  modal  density  region  Is  therefore  In  sharp 
contrast  with  the  behavior  at  high  modal  density  as  discussed  In 
Section  B. 

The  low  density  relative  variation  will  equal  the  high  density 
variation  when 

2  1  if  “V  2tn? 

%  "  (  t|  Siran;  “  ^m^  , 

or, 

2rang  -  81/32  -  2.5  .  (3.13) 

Since  2ira  Is  the  effective  bandwidth  of  a  mode,  this  repre¬ 
sents  an  overlapping  of  2  to  3  modes,  which  Is  consistent  with 
Schroeder's  criterion  for  the  passage  from  non -overlapping  to 
overlapping  behavior,  (Ref.  3).  Of  course,  the  non-overlap 
assumptions  which  have  been  made  really  preclude  our  full 
acceptance  of  any  conclusions  which  apply  the  work  of  this 
section  to  a  region  of  overlap  above  n.  .  2ra  ■  1. 

D.  STATISTICAL  ESTIMATICM  OF  THE  RESPONSE  OF  A  HARMONIC 

CSCILUTOR  ATTACHED  TO  A  FLAT  PLATE 

The  response  of  a  simple  oscillator  attached  to  a  structure 
of  many  degrees  of  freedom  will  be  approached  from  the  energy- 
statistical  method  developed  by  Lyon  and  Maldanlk  (Ref.  10).  In 
subsection  1,  we  develop  expressions  for  the  power  transfer 
betvfeen  the  oscillator  and  the  structure.  From  these  we  compute 
the  mean  emd  the  variance  of  the  power  flow,  since  as  we  Indicate 
In  subsection  2,  this  power  flow  governs  the  energy  level  of  the 
oscillator  when  It  is  attached  to  the  structure.  Using  the  known 
mean  and  variance  of  the  power  flow  In  subsection  3>  we  assume 
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an  appropriate  distribution  and  compute  confidence  levels  for  the 
vibration  levels.  It  Is  this  last  step  which  gives  engineering 
value  to  the  estimates  since  It  defines  a  "safety  factor"  In  terms 
of  so  many  db  which  must  be  added  to  the  mean  value  estimates  In 
oxvler  to  Insure  that  the  measured  levels  will  fall  below  tiie 
estimate  a  certain  fraction,  say  95^ »  of  the  time. 

1.  Power  flow  between  an  oaclllator  and  a  flat  plate 

We  assume  that  the  single  substxnicture  mode  -  many  struc¬ 
tural  mode  Interaction  problem  may  be  modelled  by  the  simple 
system  shown  In  Fig.  2.2.  The  transverse  velocity  u  of  the  flat 
plate  Is  governed  by  the  equation 


O  P_.ii 


(4.1) 


where  k  Is  the  radius  of  gyration  of  the  cross  section,  C£  Is  the 
longitudinal  plate  velocity,  (ip  Is  related  to  the  loss  factor 
i)p  for  the  plate  by  s  Pg  is  the  surface  density  of  the  panel, 

and  p  Is  the  force  per  unit  area  on  the  plate  taken  positive  In 
the  u  direction.  A  homogeneous  plate  of  thickness  h  s  2^3 
and  area  Ap  has  a  modal  density  (in  radian  freqv^ency  space) 


n  =  A^2n-c^h 


(4.2) 


The  pressure  p  v/111  be  divided  Into  the  random  noise  sources 
r  which  sustain  the  vibrations  and  the  reaction  force  from  the 
oscillator  at  the  attachment  point  Thus, 


(4.3) 


the  eigenfunctions  of  the  plate  satisfy  the  following 
relations. 


if  if  dx  »  6 


(4.4) 


29 


If  w«  expand  u,  r  and  6(2.  -  In  these  functions. 


6(x  -  5o) 


dr 

5T 


« 

m 


m 


m 


Vs) 


(^.5) 


then  Equation  (4.1)  becomes 


9 

d^J, 


dU. 


^  -ar  +  “X  +  Ps^  ’^'ra  al  “  ' 


(4.6) 


..2  _  ,,4.2.2 


Where,  a,„  -  k^x  c^  and  -  ^^(x^) 


The  oscillator  Is  shown  In  Figure  2.2  along  with  Its 
mechanical  diagram  and  mobility  analog.  If  we  consider  the  force 
on  the  mass  to  be  the  dependent  variable,  then  Its  equation  of 
motion  Is 


(4.7) 


«»>. 


Here,  1)9  -  Pq/cu  Is  the  loss  factor  for  the  oscillator,  and 
(K/M)^^  Is  the  resonance  frequency. 


30 


The  symmetry  between  Equations  (4.6)  and  (4.7)  Is  Improved 
by  the  following  change  In  variables: 


(*.8) 


Using  the  new  variables,  the  equations  become 


at“ 


(^.7a) 


du. 

dt 


“o^s(t) 


(4.8a) 


Let  us  assume  that  the  oscillator  is  excited  by  wideband 
noisej  but  the  plate  has  no  excitation  other  than  by  its  attach¬ 
ment  to  the  oscillator.  The  power  received  by  the  plate  is: 


m  m 


The  cross  moments  have  been  computed  for  the  pair  of 
Equations  (4.7a)  and  (4. to)  in  Reference  10.  Their  values  are 
(Ref.  11) 


— - 


(4.10) 


■Oie  sununatlon  (4.9)  Is  subject  to  some  uncertainties  since 
one  would  not  caz«  to  specify  or  cue  precisely  for  all  possible 
experimental  configurations.  We  ^erefore  describe  them  statletl* 
cally.  l!he  transferred  power  Is 


I? 

m 


-1 


(4.11) 


We  assume  that  the  frequencies  distributed  along  ^e 

frequency  axis  as  a  Poisson  process  with  a  mean  value  given  by 
Equation  (4.2).  The  mean  value  of  P  Is  (Ref.  8) 


n 


<»o+»p> 


r 


dm 


(*.12) 


where 

»i„  ■  *f>p'>j|i>V 


Is  the  (real)  Input  impedance  of  eui  Infinite  plate,  (Ref.  12). 
Equation  (^4.12)  may  be  considered  Idie  basis  for  mean  value  estimates 
of  power  flow  frc  a  a  shaker  Into  a  structure  and  the  vibration  of 
connected  structures  (Appendix  111,  Equation  (2.3)). 

In  order  to  find  the  variation  In  estimated  response,  we  will 
need  the  variance  of  P.  It  Is  (Ref.  9) 
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The  form  of  Equation  (4.13)  Indicates  fairly  clearly  the 
sources  of  variation  and  their  effects.  The  first  term  is  the 
effect  of  spatial  variation;  for  two  dimensional  sinusoidal 
modes  its  value  is  9/4.  The  second  factor  is  that  of  irregularity 
in  the  number  of  modes  which  the  oscillator  is  likely  to  encounter. 
The  result  is  in  the  fom  of  a  central  limit  theoremwhich  states 
that  the  standard  deviation  o  will  diminish  relative  to  the 
mean  as  the  number  of  interacting  modes 

N  «  irn(pQ+Pp)  (4.14) 


is  Increased. 

2. 


To  get  the  power  balance  equations  for  the  system,  one 
multiplies  Equation  (4.7a)  by  -JiXjn^t  and  Equation  (4.8a)  by 
-/0dt  and  averages.  The  result  is 


Relation  between 
vibration  leveis 


?ower  transfer  and  equilibrium 


4 
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In  Inference  10,  It  was  shown  that  If  0s  -  0  and  s.  Is  a  set 
of  wide  band  random  noise  sources,  then  the  power  flow  into  the 
oscillator  Is 


«rtiere  ~  ^  *  1^  we  can  assume  that  the  sources 

will  keep  all  the  structure  modes  at  the  same  energy  level 
^  -  Op  m  const,  then 


where 


or 


^  0  »  (e  -0  )p 
o  o  '  p  o'^^coup 

6  =  V  e®®"** 

coup  /  '^m 

m 


0=0  „ 
o  p  p  +E 
^  coup  *^0 


(4.16) 


where  essentially  the  sum  which  was  evaluated  as  P  In 

subsection  1.  We  note  that  If  P  «  P„,  the  statistics  of 

coup  o 

Pcoup  those  of  0Q.  In  the  event  that  »  Pq»  one  would 

have  0  .  0  and  no  variation  would  exist. 

°  P 

Our  estimates  of  variance  and  mean  values  In  the  previous 
subsection  therefore  require  equal  energy  In  the  plate  modes  If  we 
expect  then  to  predict  the  substructure  (oscillator)  response.  If 
the  Sill's  a.re  caused  by  a  point  source  at  position  then  the 
vsurtLance  In  response  will  contain  an  additional  factor  of  varla- 

tion  . 
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3.  Confidence  levels  for  estimates  of  remdom  response 

We  have  developed  expressions  for  the  mean  and  variance  of 
power  flow  or  response  when  an  oscillator  Is  attached  to  a  flat 
plate.  We  should  now  like  to  use  these  In  such  a  way  that  con¬ 
fidence  levels  for  the  estimated  response  may  be  computed.  In 
particular^  If  we  calculate  an  average  response,  then  we  might 
ask,  "by  how  many  db  should  one  Increase  this  estimate  so  that 
the  estimate  Is  not  exceeded  In,  say,  95$^  of  the  cases?"  Our 
purpose  In  this  subsection  Is  to  answer  this  question. 

We  do  not  know  the  probability  density  ♦(?),  but  If  we  did, 
the  above  requirement  would  be  met  by  setting  our  estimate  at  A, 
such  that 


0(P)  dP 


0.95  *  CL 


(^.17) 


We  shall  get  around  our  Ignorance  of  ^  by  assuming  a  form 
for  It.  We  require  a  density  function  which  Is  defined  for  posi¬ 
tive  arguments  only,  and  has  parameters  easily  related  to  the  mean 
and  variance.  These  requirements  are  met  by  the  gamma  distribution. 


„o-i  _-x/^ 

a!,  ,  (^-18) 

X  T(a) 

where  (Ref.  13) 

2  2 

m  B  Xa  ,  a  »  X  a 


One  has  then 

a  ■ 


m2 


(^.19) 


Changing  variables  to  y  ■  x/X,  the  condition  Equation  (4.17) 
becomes 
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,  rk/\  -  B 

r'  (a)  I  (Jy  y“'^  -  CL  (4.20) 

-Jo 


whArs  7(a,  B)  Is  the  Incomplete  ganms  function,  (Ref.  l4). 

The  function  Equation  (4.20)  Is  plotted  by  Jahnlce  and  Bmde 
(Ref.  15)  for  various  values  of  CL.  Using  their  values,  we  have 
plotted  curves  of  constant  CL  In  Figure  2.4,  with  the  ratio  of 
squared  mean  to  variance  as  the  abscissa,  and  10  log  A/m  as  the 
oz^Unate,  where 

k/a  mB/a  .  (4.21) 

The  curves  In  Figure  2.4  essentially  tell  us  how  many  db  to 
add  to  the  BMan  value  estimate  In  order  to  ensure  that  a  fraction 
CL  of  the  measured  values  of  response  will  fall  below  the  estimate. 
Thus  for  example.  If 

N  -  rn  (pQ+Pp)  -  10  , 

then  a  -  4.44  and  one  would  have  to  add  2.6  db  for  a  confidence 
level  of  95$(  or  3*8  db  for  a  confidence  level  of  99ft,  If  the 
plate  Is  driven  by  a  point  source  shalcer  so  that  an  additional 
variance  factor  arises  from  the  modal  energy  uncertainty,  ttien 
the  relevant  value  of  a  is 


Referring  to  Figure  2.4,  a  confidence  level  of  9^  would 
require  that  we  add  3*9  <lb  to  the  mean  value  estlmte.  For 
confidence,  we  would  have  to  add  5*2  <lb* 
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FIG.  2.2  DIAGRAM  OF  PLATE  WITH  ATTACHED 
OSCILLATOR 
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SHAKER  FORCE 
f,(t) 


FIG.  2.3  OSCILLATOR  SKETCH,  MECHANICAL  DIAGRAM, 
AND  MOBILITY  ANALOG  EQUIVALENT  CIRCUIT 
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10  LOG  A/m 
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FI6.2.4  SAFETY  FACTOR  IN  ESTIMATION  REQUIRED  TO 
ACHIEVE  GIVEN  CONFIDENCE  LEVEL 
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APPENDIX  Ill 


ENERGY  SHARING  BETWEEN  CONNECTED  STRUCTURES 


A.  INTRODUCTION 

When  large,  complex  structuz*e8  are  deslmed  to  operate  In 
environments  that  may  cause  them  to  vibrate  (e.g..  Intense  sound 
fields,  boundary  layer  turbulence,  etc.)^  design  engineer 
Is  usually  concerned  with  the  response  of  a  particular  substruc¬ 
ture  that  may  be  susceptible  to  malfunction  (e.g.,  components  of 
the  structures  that  carry  electronic  equipment) .  We  wish  to 
develop  methods  for  estimating  the  vibration  of  substructuz^s 
attached  to  large  complex  structures.  We  define  a  large  complex 
structure  such  that  It  has  high  modal  density,  and  Its  modes  are 
fairly  uniformly  distributed  In  frequency  space.  We  assume  that 
the  substructure  has  Its  own  modal  density  and  that  the  method  of 
attachment  of  the  substructure  to  the  large  structure  Is  known. 

We  further  assume  that  the  vibration  of  the  structure  Is  pre¬ 
dominantly  resonant,  l.e.,  the  modes  of  vibration  of  both  the 
structure  and  substructure  are  not  critically  damped.  Under  this 
assumption  one  may  treat  the  Interaction  between  the  structures 
by  considering  the  Interaction  of  one  set  of  oscillators  (modes) 
that  Is  coupled  In  a  known  manner  to  another  set. 

The  power  flow  between  two  linearly  coupled  oscillators  has 
been  studied  by  Lyon  and  Maldanlk  (Ref.  1) .  The  formalism  has 
been  extended  to  cover  the  case  of  the  power  flow  between  two  sets 
of  coupled  oscillators.  The  theory  has  also  been  applied  specifi¬ 
cally  to  the  Interaction  between  a  reverberant  acoustical  field 
(considered  to  have  a  very  high  modal  density),  and  a  structure, 
e.g.,  a  beam,  a  panel  or  a  combination  of  both,  (Ref.  2).  It  Is 
possible  to  treat  the  connected  structu.’e  problem  by  the  same 
methods,  as  we  shall  show. 

The  coupling  Is  strongest  between  modes  of  the  structure  (high 
modal  density  assumed)  and  a  mode  of  the  substructure  when  they  lie 
close  In  frequency  space.  Modes  which  lie  outside  each  others  pass- 
bands  are  usually  weakly  coupled.  In  a  given  frequency  band 
centered  about  a  given  a>  (to  =  angjular  frequency).  If  ©as  the 
modal  energy  of  the  substructure  and  9^  Is  the  modal  energy  of  the 
structure,  then  If  the  structure  la  excited  by  external  source 
fields,  the  relation  between  Sgg  and  63  can  be  shown  to  be  (Ref.  1) 


ss 
s  ''s 


9.  p. 


(1.1) 


where 
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(1.2) 


.88  'C 

”  V%8 

The  ratio  (1.2)  la  called  the  coupling  factor  for  the  sub¬ 
structure;  and  T)ss  loss  factors  associated  with  the 

coupling  and  the  substructure  damping  respectively.  While  1)55 
Is  related  to  the  power  dissipated  Internally,  tIq  Is  related  to 
the  power  lost  (or  gained)  by  the  substructure  to  (or  from)  the 
structure  via  the  coupling  (c.f.,  R__j  =  aM_  .  for  the 
acoustic  case  (Ref.  l)).  P 

Equation  (1.1)  states  that  In  the  steady  state  condition  the 
modes  of  the  substructure  attain  modal  energies  that  are  a  fraction 
a  of  the  modal  energy  of  the  structure.  For  ®- 

and  for  Tjc«n„„,  e„«  e„.  ®  ® 

Equation  (1.1)  Is  an  approximate  statistical  formula,  (Ref.  1) . 
The  approximation  Is  better  as  the  number  of  Interacting  modes  In 
a  given  frequency  band  Increases.  This  Is  so  because  the  formula 
Is  derived  by  averaging  a  summation  over  modes  that  lie  In  the 
frequency  band  (e.g.,  octave  or  third-octave) .  Thus  one  would 
expect  the  approximation  to  be  better  the  higher  the  modal  density 
of  the  structures  and  the  wider  the  frequency  band  chosen.  One 
could  get  an  Idea  of  the  errors  involved  In  the  smoothing  pro¬ 
cedures  by  computing  the  variances  of  the  fluctuations  In  the 
responses  of  the  structure  and  substructures;  this  Is  being 
reported  In  Appendix  11.  In  this  Appendix  for  the  present  we 
concentrate  on  examining  whether  the  general  assumptions  and  Ideas 
expressed  with  respect  to  Equation  (1.1)  are  In  keeping  with  exper¬ 
imental  data. 

B.  BEAM  CANTILEVERED  TO  FLAT  PLATE 

To  test  Equation  (1.1),  we  choose  to  examine  a  beam  canti¬ 
levered  to  a  flat  plate.  The  beam  and  the  plate  are  so  chosen 
that  the  coupling  loss  factor  t}q  may  be  conveniently  estimated. 
Moreover,  the  dimensions  of  the  beam  and  the  plate  are  such  that 
the  conditions  for  the  validity  of  Equation  (1.1)  are  approxi¬ 
mately  satisfied.  Below  we  discuss  briefly  the  theoretical  and 
experimental  results  Involved  In  the  test. 

1.  Theoretical  Considerations 


To  determine  t)q,  we  assume  that  the  beam  Is  excited  by  an 
external  source,  and  we  calculate  the  power  loss  to  the  plate 
via  the  coupling.  The  power  flow  Into  the  plate  from  the  beam 
Is  given  by 
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M  Re 


» 


(2.1) 


where  Zp  Is  the  moment  Impedance  of  the 

beam  at  the  Junction  and  is  the  mean 
at  the  Junction  by  the  vibration  of  the 


plate  as  seen  by  the 

square  moment  generated 
beam. 


In  terms  of  n  ,  the  power  that  the  beam  loses  to  the  plate 
c 

through  the  coupling  is  given  by 


^p  =  \  “  «b 


(2.2) 


where  is  the  mean  square  velocity  of  the  beam  averaged  in  time 
and  space  and  is  the  total  mass  of  the  beam. 

Equating  Equations  (2.1)  and  (2.2)  we  have 


% 


_z_ 


(2.3) 


If  a  typical  dimension  of  the  Junction  is  small  compared  to 
a  wavelength  in  the  beam  and  plate,  then  (Ref.  2) 


Re 


P  P  P 
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where  pg  Is  the  density,  hn  Is  the* thickness, 
longitudinal  velocity  of  the  plate. 


(2.4) 


and  Cp  Is  the 


,  We  have  now  to  determine  M  .  Consider  the  analog  circuit 
of  Figure  3.1  depicting  the  coupling  between  the  plate  and  beam. 
The  flow  quantity  Is  the  angular  velocity  and  the  potential  drop 
Is  the  moment  M.-ZM  and  are  the  moment  Impedances  of  the 

plate  and  the  beam  as  seen  by  the  beam  and  plate  respectively. 

We  note  from  Figure  3.1  that  the  mean  square  moment  Is  given  by 
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Zw  +  Z‘‘ 

D  P 


(2.5) 


-gb  -gb 

Where  M  Is  the  mean  square  moment  of  the  beam,  M  is  related 

to  the  mean  square  velocity  on  the  beam,  ^  by  the  equation 


(2.6) 


where  pw  is  the  density,  Cb  the  longitudinal  velocity,  the 
radius  of  gyration  and  Sb  the  cross-sectional  area  of  the  beeun. 

Combining  Equations  (2.3),  (2.4),  (2.5),  and  (2.6)  we  obtain 


n  =  -=V- 

°  4h^ 


(2.7) 


where  hb  Is  the  thicloness,  w  the  width,  i  the  length  of  the  beam, 
emd 


b  p 


The  expression  for  4  may  be  written  In  the  form  (Ref.  2) 


where 


Pb^b®b^''S) 

32p,h3;f  ' 

P  P  P 
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Ic^  being  the  vraivenumber  for  the  beaun.  For  a  olTOular  beam  of 

radius  r,  r  ■  log^  (kpr)  ikpr«l],(Ref .  2).  For  a  rectangular 

beaun  the  expression  for  T  Is  lacking.  We  assume  that  for  our 
arraungement  (see  Figure  3*2)  the  absolute  value  of  P  Is  t2ie  order 
of  unity.  Thus,  for  the  case  under  consideration  here 

®p"®b^  "  “  ^  cm  amd  k^«>0.7  cm~^  at  4  kc),  we  may  set  i-1 

In  the  frequency  range  from  50  cps  to  12  kc.  One  should  note  that 
a  better  approximation  may  lead  to  somewhat  lower  values  of  ^  and 
hence  smaller 

2.  Experimental  Considerations  and  Results 

Equation  (l.l)  as  It  Is  given  Is  not  convenient  for  practical 
purposes  since  one  cannot  measure  dgg  and  dg  directly.  We  have  to 
express  the  moda.i  ener^es  In  terns  of  measurable  qusuitltles.  The 
energy  of  a  structure  (^a  beam  or  a  plate)  In  a  given  fraquency 
band  (e.g.,  third-octave  or  octave)  Is  given  by 


where  S^ls  the  acceleration  spectral  uenslty  of  the  structure 
In  the  frequency  band  /ku,  <a  is  the  center  frequency  of  the  band 
and  Ms  Is  the  total  mass  of  the  stznicture.  Ihe  modal  energy  of 
the  vibration  of  ttie  structure  Is  given  by 


S  ^s  MS^ 


(2.9) 


where  n(ca)  Is  the  average  modal  density  In  the  band  da>. 

Substituting  Equation  (2.9)  In  Eqviation  (l.l)  and  making 
the  appropriate  Identification  of  the  quantities,  we  obtain 


ga  ”  np(a))M^  \  +  %  ' 


(2.10) 


where  subscripts  p  and  b  refer  to  quantities  related  to  the 
plate  and  the  beam  respectively. 

Equation  (2.10)  Is  now  In  a  form  where  the  parameters  Involved 
are  measurable. 
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The  beam-plate  system  is  shown  In  schematic  form  In  Figure 
3.2.  The  dimensions  of  the  aluminum  beam  are  95"  x  1.5”  x  1/8" 
and  those  of  the  plate  72"  x  48"  x  1/8".  The  choice  of  a  rather 
thick  beam  and  plate  Is  made  to  avoid  excessive  loading  of  the 
structures  by  the  accelerometers  at  high  frequencies.  The  beam 
Is  attached  perpendicularly  to  the  plate  by  epoxy.  The  positions 
on  the  plate  and  beam  where  measurements  were  taken  Is  Indicated 
In  Figure  3.2  by  the  circled  numbers. 


llie  modal  density  of  the  plate,  np.  Is  measured  by  exciting 
the  plate  with  a  shaker  and  counting  pure  tone  resonances.  It 
Is  found  that  the  modal  density  of  the  plate  conforms  closely  to 
Its  asymptotic  theoretical  value,  which  Is  given  by  (Ref.  3) 


np(a)) 


VT  Ap 

arc  V 
p  p 


(2.11) 


where  Ap  Is  the  area  of  the  plate. 

For  our  plate  np(a>V2^3.8  x  10  modes/radian  per  second. 
Deviation  from  the  value  of  np(a>)  as  given  by  Equation  (2.11) 
occurs  only  below  about  30  cps.  We  shall  not  be  concerned  with 
frequencies  below  50  cps. 


The  modal  density  of  the  beam  is  measured  in  the  same  way. 
It  is  found  that  the  measured  modal  density  of  the  beam  agrees 
fairly  well  with  the  asymptotic  theoretical  value  given  by  the 
formula  (Ref.  3) 


np(<o) 


i  (V3 

1“^  I 


(2.12) 


It  is  of  interest  to  note  that  above  6  Vcc  the  modal  density 
of  the  beatm  exceeds  the  value  which  is  given  by  Equation  (2.12) 
and  slowly  approaches  the  modal  density  of  a  plate  as  the  frequency 
is  increased.  From  these  measurements  the  quantity 


np(a))M 
^°®10  npI<o)M 


£ 

b 


(2.13) 


Is  determined.  This  quantity  Is  given  In  graphical  foimi  In  Figure 
3.4.  The  factor  given  by  Equation  (2.13)  represents  the  ratio 
between  the  acceleration  spectral  densities  of  the  plate  and  the 
beam  (in  db)  for  the  case  of  "thermal  equilibrium"  l.e.,  equl- 
partltlon  of  modal  energy  (0._  ■  0.). 

BS  8 
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llie  quantity  %  Is  determined  by  measuring  the  reverberation 
time  of  the  beam  (unattached) (Ref .  4).  It  Is  related  to  the 
reverberation  time  T]g  by  the  equation 


The  loss  factor  Is  determined  for  two  cases,  an  undamped 
beam  and  a  damped  beam.  The  damping  was  achieved  by  attaching 
damping  tape  to  the  beam.  Using  the  theoretical  values  of 
[Equation  i2.7)  with  4=1]  and  the  experimental  values  of 
^uatlons  (2.13)  and  (2.14),  S^S|  Is  determined  from  Equation 

(2.10) .  The  results  are  given  In  graphical  form  In  Figures  3.3 
and  3.4. 

We  then  proceed  to  measure  the  quantity  s|/Sp  directly  In 
order  to  compare  It  with  the  values  obtained  above.  This  con¬ 
stitutes  a  test  of  the  validity  of  Equation  (2.10). 

A  block  diagram  of  the  system  for  measuring  sf/Sn  Is  shown 
In  Figure  3.5.  The  acceleration  Is  measured  (in  third-octave  and 
octave  bands)  successively  at  selected  points  on  the  plate  and 
beam  (see  Figure  3*2)  and  simultaneously  at  point  P  on  the  plate. 
The  difference  (in  db)  between  the  acceleration  at  a  given  point 
and  at  point  P  Is  thus  obtained.  In  a  given  frequency  band  these 
differences  are  averaged  over  all  the  selected  points  of  the  beam 
and,  separately,  over  all  the  selected  points  on  the  plate.  The 
difference  between  these  averages  gives  the  value  of  In  the 

frequency  band.  The  values  of  S|/Sp  so  obtained  are  given  In 
graphical  form  In  Figures  3.3  and  3.4  for  two  cases  of  damping; 
the  two  cases  are  the  same  as  those  for  which  reverberation  time 
measurements  were  made. 

It  should  be  pointed  out  that  the  beam  and  the  plate  are 
found  to  be  fairly  3?everberant  In  the  sense  that  the  accelerations 
(normalized  by  the  acceleration  at  point  P)  taken  at  the  selected 
points  on  the  plate  and  the  beaun  are  found  to  be  within  1  or  2  db 
of  their  averaged  values.  Note  that  this  Indicates  that  the  vibra¬ 
tion  of  the  structures  Is  predominantly  resonant,  a  necessary 
condition  for  the  validity  of  Equation  (2.10). 

The  cunparlson  between  the  direct  measurement  of  st/Sp  and 
the  Indirect  measurement  (where  the  right-hand  side  of  Equation 

(2.10)  Is  computed)  Is  shown  graphically  In  Figures  3.3*  and  3.4. 

We  feel  the  agreement  between  the  values  v/ould  be  satisfactory 
as  an  engineering  estimate.  The  point-to-point  variations  are 

In  some  cases,  however,  quite  large.  It  was  pointed  out  earlier 
that  when  the  number  of  modes  either  In  the  plate  and/or  the  beam 
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that  participate  In  the  coupling  In  any  given  frequency  band  Is 
not  large,  fluctuations  from  the  Ideal  case,  which  Equation  (2.10) 
represents,  are  to  be  expected.  Although  the  number  of  modes  In 
the  plate  In  any  given  third-octave  band  Is  fairly  large,  the 
niunber  of  beam  modes  Is  rather  small.  It  Is  Interesting  to  note 
that  the  fluctuations  are  greatly  suppressed  when  measurements 
are  conducted  In  octave  Instead  of  third-octave  bands. 

The  results  obtained  Indicate  that  the  method  developed  above 
for  estimating  the  vibrations  of  substructures  may  prove  fruitful. 
Some  refinements,  especially  with  respect  to  estimating  variances 
from  the  average  values,  would  prove  useful. 

C.  THE  VIBRATION  OF  CONNECTED  PANELS 

1.  Introduction 


We  shall  now  conslaer  a  connected  structure  which  Is  a 
typical  part  of  electronic  paclcages.  The  racks,  housings,  and 
chassis  of  electronics  systems,  for  exstmple,  are  made  up  of 
plates  or  peuiels.  We  ask  how  much  a  subpanel  will  vibrate  when 
Its  mounting  rack  Is  randomly  excited.  In  order  to  answer  this 
question  eventually,  we  Investigate  here  the  flow  of  vibrational 
power  between  plates  which  are  attached  to  each  other  by  welded 
seams.  We  do  not  Intend  to  e^qplore  each  of  the  great  number  of 
Joints  between  plates  which  are  In  use  but  rather  to  apply  the 
general  methods  which  have  been  previously  outlined  to  a  few 
types  of  connection. 

Our  fundamental  fomtula  Is  the  energy  ratio  of  two  coupled 
modes,  given  In  Equation  (l.l).  If  the  structures  are  reverberant, 
as  we  shall  assume,  we  can  eliminate  the  modal  energies  In  favor 
of  the  "energy  velocity"  spectral  densities  and  obtain 


"ss 


(3.1) 


where  n  Is  the  modal  density  smd  M  the  mass  of  the  respective 
structure. 


For  a  flat  plate,  the  modal  density  Is 


(3.2) 
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where  A  is  the  (one-side)  surface  area,  h  the  thickness  of  the 
plate,  and  c^  Is  the  plate  longitudinal  wave  velocity.  The  mass 
of  the  plate  Is  M  »  pAh,  where  p  Is  the  mass  density.  Hence, 
one  obtains  for  connected  plates 


ss 


ph^c^ 

8 

2  1 
ph  c^ 

|ss 

ss 


(3.3) 


For  plates  of  the  same  material  we  have  simply 


ss 


(3.4) 


This  formula  allows  us  to  calculate  the  vibration  level  of 
the  substructure  If  we  knov;  the  thickness  ratio  and  the  coefficient 


SS 

The  deteinnlnatlon  of  p,  will  occupy  us  In  the  next  two 
sections.  If  the  substructure  Is  excited,  then  pis  may  be  defined 
In  terms  of  the  power  It  loses  by  transmission  and  dissipation  In 
the  following  manner 


ss  _  Power  transferred  to  structure 
^s  'Total  power  lost 


(3.5) 


If  these  quantities  are  Interpreted  In  terms  of  loss  factors. 
Equation  (l.l)  results.  Besides  other  results,  we  will  present  In 
subsection  2.  theoretical  expressions  for  the  coupling  loss  factor 
of  a  plate  which  Is  clamped  to  a  second  plate.  In  subsection  3. 
we  will  report  related  measurements.  The  technique  used  In  sub¬ 
section  2.  Is  wave  theory.  The  plates  are  taken  as  Infinitely  large 
and  a  reverberant  field  Is  simulated  by  averaging  over  all  direc¬ 
tions.  Ihe  boundary  absorption  equations  used  In  subsection  3. 
are  derived  from  a  plate  analogy  to  room  acoustics.  The  absorption 
Is  obtained  by  measuring  the  reverberation  time. 


51 


2.  Flexural  Wave  Tremamlsslon  between  Clamped  Panels 

We  consider  a  thin  Infinite  half -plate  1  Joined  under  a 
right  auigle  to  a  thin  Infinite  plate  2  (Figure  3*6).  We  assume 
that  the  plates  are  Joined  as  If  they  vere  made  from  one  piece. 
The  Joint  may  be  approximated  by  a  continuous  welded  or  solde3?ed 
seam. 


Let  u,  V,  v,  be  the  displacements  In  the  directions  x,  y,  z. 
For  small  displacements  v  ■  0.  At  the  Joint  x  =  0,  z  »  0,  the 
transverse  displacements  of  both  plates  must  vanish,  u  >  0, 
w  B  0,  (u  refers  to  plate  2  amd  w  to  plate  1,  but  It  Is  not 
necessary  to  use  Indices) .  The  conservation  of  the  right  angle 
demands  that  dw/dx  >  -du^z.  Let  Q  be  the  moment  (per  unit 
length)  acting  on  plate  1  at  the  Joint  (Fl^re  3.7;. 

From  Love  (Ref.  5)  the  bending  moment  along  the  straight 
edge  of  a  plate  Is 

G  =  -DLS^w/Sn^+od^w/Ss^J  (3.6) 

where  n  Is  the  coordinate  out  of  the  edge  (s  x)  and  s  the 
coordinate  along  the  edge(s  y),  and 

D  B  pc^Vl2  (3.7) 


Is  the  flexural  stiffness  of  the  plate.  Applying  this  to  each 
of  the  three  half  plates  we  obtain 


=  D2(d^u/5z2)  -  D2(d^u/a2^)+0  (3.8) 


Our  set  of  boundary  conditions  Is  now  complete. 

The  differential  equation  for  small  transverse  vibrations 
of  a  thin  plate  Is  [DV^  +  phd2/5t2]  w  b  o.  For  harmonic  time 
dependence  exp  lot  suid  exp( -Ikyslna)  on  y  we  get  the  solutions 

exp(+lkxco8a)  and  exp(+kx  Vl+sln^)  In  x.  We  calculate  the 
traninilsslon  of  a  single  wave 

w  B  exp(la)t  -  ik^yslna  -  Ik^xcosa)  ,  (3.9) 
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Incident  in  plate  1  at  an  angle  of  incidence  a  (Figure  3.6). 
Satisfying  the  boundary  conditions  we  find  for  the  wave  trans¬ 
mitted  into  plate  2, 


u  »  C  exp(iOi)t  -  ikC^ysinP  -  Ik^zcosP)  (3.10) 

where 

C  *  -(1+A)  k^r/kg  (3.11) 

»  -2i(k^r/k2)cosa[icosa  -  Vl+sin^  +  r(icosp  -Vi+sin^p) 

r  *  DQ^k^/2D2k2  >  (3.12) 

and  A  is  the  amplitude  of  the  reflected  wave.  To  preserve  the 
dependence  on  y  the  angles  a  and  p  are  related  by 


k^sina  s  kgSinp  (3.13) 

The  mean  square  velocity  on  plate  2  is 


2  2 

(a>k^r/k2)  cos  a 
l+rVi  -sin^aVl  -(  k^k2)  ^sii  j^a+r  A+3Tn^(^l-i-(  k^/k2y^sin^a+r^ 


Vr.  * 


0) 


(3.14) 

If  the  transmission  from  plate  1  into  plate  2  is  small,  we 
can  take  the  amplitude  of  the  reflected  wave  as  A  =  -1,  so  that 
the  squared  velocity,  averaged  over  plate  1,  is 


(3.15) 


which  is  independent  of  a. 
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In  a  reverberant  and  diffuse  field  waves  come  from  all 
directions  with  equal  probability.  The  angular  average  of 
Equation  (3>1^}  oa.n  be  readily  found  In  two  special  cases.  For 
equal  plates 


<v|>(j 


1 


(3.16) 


and  for  plate  1  much  thinner  than  plate  2  (k^»lC2j  Dj^«D2,  r«l). 


<4>a 


1 

5F 


^2  ^2 
2arcsin  ^  +  sln2arcslrj  ^ 


*^1  ^1 

2rk|D| 


(3.17) 


The  average  power  per  unit  boundary  length  carried  by  a  wave 
incident  at  the  angle  a  (with  amplitude  1)  equals  the  average 
energy  density 


E  -  2pha)^  (3.18) 

times  the  group  velocity 

Cg  -  [dk/doj]’^  »  2u)/lc  (3.19) 

times  cosa  .  The  corresponding  power  transmitted  to  plate  2 
(and  hence  lost  to  plate  1)  Is 


^’lost  “ 


(3.20) 


The  ratio  of  the  angular  averages  is  the  absorption  coefficient 
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r»/2 


coB^a  do 


j  •>r/2  l+r -8ln^o''4  -( k^^/kg)  ^Bin^^r"^  l+Bln^ayi+(  k^Ag)  ^Bin^^r* 


(3.21) 

which  can  be  calculated  readily  in  two  special  cases.  (The 
integrand  is  to  be  taken  as  zero  where  it  is  complex.)  For 
equal  plates  we  find 


•y  =  V27 


(3.22) 


and  for  plate  1  much  weaker  than  plate  2 


_  I"  kg  kg 

7  =  IT-!  2  arcsln  i—  +  sln2arcsln 

*  ^  L  1  1 


(3.23) 


{  de  ai^e  looking  at  plate  2  as  a  boundaxT-  absorber  attached 
to  plate  1.  Note  that  this  step  is  an  Intermediate  step  in  the 
calculation  of  which  In  turn  Is  needed  to  calculate  the  excita¬ 
tion  of  plate  1  when  plate  2  la  driven^ 


It  remains  to  relate  7  to  t)c«  This  Is  easily  done  using  the 
formula,  (Ref.  6)  for  the  reverberation  time  T 


2.2 

ir 


yc  h 


(3.24) 


where  L  la  the  length  of  the  absorbing  boundary.  We  may  inter¬ 
pret  the  transmission  as  a  coupling  loss  factor 


%1 


oj  A, 


(3.25) 


which,  together  with  the  Internal  loss  factor  t)i  allows  us  to 
compute  the  coupling  factor 
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(3.26) 


with  plate  1  assuming  the  role  of  the  substructure  and  plate  2 
the  role  of  the  structure. 

3.  Experimental  Study 

Measurements  were  made  on  one  aluminum  plate  and  on  two 
steel  plates  ("half  plates")  Joined  to  heavier  plates  of  similar 
material.  The  light  plates  were  cut  Irregularly  to  promote 
diffuse  conditions.  The  heavier  plates  were  damped  to  approxi¬ 
mate  Infinite  plates  when  seen  from  the  attached  light  plate. 

The  light  plate  was  struck  and  the  resulting  vioratlon  was 
sensed  with  an  accelerometer  (Glennlte  31^).  After  filtering 
through  a  third-octave  band  filter  set  (B+K  I609)  the  decay  time 
was  measured  (Ref.  4)  (bbN  Reverberation  Time  Meter).  The 
recorded  times  were  converted  Into  loss  factors  and  absorption 
coefficients  according  to  Equations  (3.25)  and  (3.26).  The  results 
are  plotted  In  Figures  3.8-10  together  with  theoretical  values  of 
y  from  Equation  (3.26). 

We  started  with  a  thin  aluminum  plate  which  we  glued  to  a  heavy 
one  by  means  of  an  artificial  resin  (Fig.  3.8).  Tlie  measured 
values  of  the  absorption  coefficient  turned  out  to  be  large  In  com¬ 
parison  with  theory.  We  suspected  that  some  Imperfections  In  the 
Joint  or  the  glue  Itself  was  causing  the  excess  absorption.  This 
suspicion  was  confirmed  when  we  measured  the  total  loss  factor  q' 
of  the  same  plate  glued  to  a  heavy  steel  bar.  The  difference 
between  q'  and  the  free  plate  loss  factor  q,  represents  the  losses 
of  the  Joint  and  Is  considerable. 

We  consequently  went  from  the  glued  Joint  to  an  all-metal 
Joint.  Because  of  the  large  size  of  the  plates  we  had  difficulties 
In  the  production  of  the  test  object.  The  next  plate  we  tested 
was  a  thin  steel  plate  welded  to  a  heavy  one  (Figure  3.9). 

The  welded  Joint  was  not  uniform  and  we  found  again  absorption 
values  which  were  higher  than  the  predicted  ones,  v/e  finally 
obtained  a  very  good  weld  on  a  steel  plate  Joined  to  another  one 
about  three  times  as  thick.  In  this  case  the  agreement  between 
experiment  and  theory  was  excellent.  It  seems  that  the  slightly 
Increased  torsional  stiffness  of  the  Joint  (due  to  the  welded 
depositions)  did  not  make  the  theory  unappllcable . 
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Oood  agz>eement  between  the  presented  theory  auid  experiments 
cannot  be  expected  at  low  frequencies  where  the  Joint  Is  only  a 
few  flexural  wavelengths  long  emd  also  at  high  frequencies  (above 
coincidence)  where  radiation  dominates  the  Internal  loss  factor. 


Vibratory  Response  of  Cantilevered  Panel 


In  subsection  1,  we  derived  an  expression  shown  In  Equation 
for  the  vibratory  response  of  a  panel  substructure  attached 
to  another  randomly  excited  peuiel  structure.  We  should  like  to 
apply  this  formula  to  the  welded  plate  configuration  shown  In 
Figure  3*10.  In  terms  of  the  thicknesses  h^  and  h2  of  plates  1 
and  2  respectively  and  the  coupling  and  loss  factors,  the  acceler- 
tion  spectrum  ratio  Is 


For  the  plates  studied,  h^  =  .057"  and  h2  *  .118".  The 
coupling  He  la  determined  from  the  measured  value  of  y  given  In 
Figure  3.10  by  using  Equation  (3.26).  The  substructure  damping 
T)i  was  measured  v;lth  and  without  applied  damping.  A  comparison 
of  the  measured  response  ratio  with  the  computed  values  without 
applied  dsunuplng  on  plate  1  is  shown  In  Figure  3*13.  A  similar 
comparison  Is  shovm  v/lth  plate  1  (substructure  panel)  more  heavily 
damped. 

Although  the  agreement  achieved  Is  not  precise,  the  values 
tend  to  lie  within  a  3  db  range  of  each  other  and  we  would  pjrabably 
regard  the  calculations  as  a  reasonable  engineering  estimate  of 
the  observed  response. 

D.  REDUCTION  OF  BEAM-PLATE  COUPLING  BY  NOTCHING  BEAM 

The  beam  which  was  described  In  Figure  3.2  was  notched  for 
the  purpose  of  reducing  the  plate -beam  coupling.  The  location  of 
the  notch  and  definitions  of  dimensions  are  shown  In  Figure  3.11. 
Since  the  reduction  of  coupling  Is  only  effective  In  reducing  the 
vlbx*atlon  when  tje  <  qjj,  the  experiments  were  run  with  applied 
damping  on  the  beam.  The  octave  band  results  with  no  notch  were 
shown  In  Figure  3.^  and  are  reproduced  In  Figui*e  3.12  under  the 
caption  w  ■  1.5"»  ■  0. 
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The  width  of  the  Joining  stud  Is  Wg  and  Its  length  Is  ig. 

The  ratio  of  beam  to  plate  vibration  Is  shown  for  various  values 
of  these  dimensions  In  Figure  3*12.  Qualitatively  the  longer  and 
thinner  the  stud,  the  greater  is  the  reduction  In  vibration.  We 
do  notice  a  resonance  effect  for  the  stud  ig  ■=  0.8",  Wg  =  O.l" 
at  about  6  kc.  This  corresponds  to  k^s  ^^1  where  k  Is  the  wavenumber 
for  flexural  waves. 
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FIG.3.1  ANALOG  CIRCUIT  FOR  THE  MOMENT 
AT  THE  BEAM-PLATE  INTERFACE 
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FIG.  3.2  SCHEMATIC  DIAGRAM  OF  BEAM-PLATE  SYSTEM 
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THIRD  OCTAVE  BAND  CENTER  FREQUENCY  IN  CYCLES  PER  SECOND 

FIG. 3.3  DWECT  AND  INDIRECT  MEASUREMENTS  OF  RESPONSE  RATIOS  WITHOUT 
APPLIED  DAMPING 
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FIG.  3.5  EXPERIMENTAL  FLOW  DIAGRAM  FOR  DIRECT 
MEASUREMENT  OF  S?/S® 
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FIG.  3.6  PLATE  GEOMETRY 
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FIG.  3.7  PLATE  BOUNDARY  CONDITIONS 
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THIRD  OCTAVE  BAND  CENTER  FREQUENCY  IN  CYCLES  PER  SECOND 
FIG.  3.9  LOSS  FACTORS  AND  ABSORPTION  COEFFICIENT  OF  0.031“  WELDED  STEEL  PLATE 


FIG.  3 
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FIG. 


70 


FIG.  3.  12  THE  RESPONSE  OF  A  BEAM  CANTILEVERED  TO  FLAT  PLATE 
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-o-ost”  thickness 
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(PLATE  1  LIGHTLY  DAMPED) 


